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A new mixing layer can be generated if the rotation of either of the two cylinders in a Taylor-Couette apparatus varies 
discontinuously along the symmetry axis. The mixing zone between the two resulting co-current streams gives rise to 
radial vorticity in addition to the primary axial vorticity. An analytic solution for the azimuthal velocity has been 
derived from which we show that the width of the mixing zone varies only with radial position. 
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1. Introduction 

A Taylor-Couette (TC) apparatus is a popular device 
in which the flow of a viscous fluid between two concen- 
tric circular cylinders is driven solely by the rotation of 
one or both of the cylinders. The TC configuration was 
probably first used to determine the viscosity of water by 
measuring the torque on the fixed cylinder [TJ [5]. TC- 
apparati are still in daily use in viscometry. According to 
his expectations Couette found that the measured torque 
increased linearly with the angular velocity at low rota- 
tion rates [2] • The deviation observed at somewhat higher 
rotation rates was ascribed to a transition to turbulence. 
Taylor [3J, however, suggested that the unidirectional Cou- 
ette flow became linearly unstable and an array of toroidal 
vortices evolved. The existence of such pairs of counter- 
rotating vortices, now known as Taylor vortices, was also 
verified experimentally by Taylor [3j. At even higher ro- 
tation rates, a wealth of different flow regimes may occur, 
see e.g. Refs. |H [5]. 

The geometrical simplicity of the TC configuration in 
combination with the rich physics of the resulting flow 
field make it popular among applied mathematicians and 
physicists to explore hydrodynamic stability, transition, 
and turbulence. A TC-apparatus offers an abundance of 
flow features in a small and closed system which makes 
it particularly appealing to experimentalists. The experi- 
mental investigation of migration of spherical particles sus- 
pended in a Newtonian fluid by Tetlow et al. [6] may serve 
as an example of a particular practical application. 

In a recent study, Sprague et al. [7] manipulated the 
inner cylinder in a TC-apparatus with a view to tailor- 
make the Taylor vortices. Two different concepts were 
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explored by means of spectral-element simulations. In the 
type A configuration the radius of the inner cylinder var- 
ied discontinuously whereas in the type B configuration 
the angular velocity of the inner cylinder exhibited a dis- 
continuous change. In this way adjacent regions of stable 
and unstable flow (with respect to the formation of Taylor 
vortices) could occur, or Taylor vortices of different wave- 
length on the two sides of the discontinuity could form. 
A related experiment with continuously changing cylinder 
profiles followed [8]. 

In the present paper we are concerned only with the vis- 
cous flow which arises at relatively low rotation rates, i.e. 
before any toroidal Taylor vortices develop. Nevertheless, 
the flow will be referred to as a Taylor-Couette flow in 
order to distinguish it from the plane Couette flow. We 
consider a configuration in which the rotating cylinder is 
divided in two parts which rotate steadily with different 
angular velocities. The differently rotating parts of the 
cylinder give rise to adjacent flow regions with different 
azimuthal velocities. An explicit analytic solution will be 
derived for the steady flow in the mixing zone between the 
two co-current streams and the axial width of the mixing 
zone will be determined. Comparisons will also be made 
with a recent analytic solution for two co-current plane 
Couette flows [9]. 

2. Mathematical formulation 

Let us consider the incompressible flow of a Newtonian 
fluid. The fluid motion in a TC-configuration is conve- 
niently described in a cylindrical coordinate system (r, 9, z) 
and V r , Vg, and V z are the velocity components in the ra- 
dial, azimuthal, and axial directions, respectively. If the 
fluid motion is sufficiently slow, and in the absence of end- 
wall effects, the velocity components in the radial and ax- 
ial directions are assumed to be negligible. We are only 
concerned with steady-state flow. Due to the axial sym- 
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metry, the unidirectional fluid motion V$ is independent 
of 9. The governing equation of motion in the azimuthal 
direction thereby simplifies to: 

\rd r + d 2 z - AW = 0. (f ) 

r r z J 

Even though the velocity component in the radial direction 
V r is assumed to be negligibly small, the radial component 
of the Navier-Stokes equations becomes: 

-V e 2 - -d r P = (2) 
r p 

where p is the fluid density. This is simply a balance be- 
tween the pressure force set up in the radial direction and 
the centrifugal force associated the azimuthal fluid mo- 
tion. The pressure field does not affect the fluid motion 
and P can therefore be obtained from Eq. (§ as soon as 
the azimuthal velocity field V$ has been determined by 
integration of Eq. ([lj . 

The fluid motion in the annular gap between inner and 
outer cylinders with radii a and b, respectively, is driven 
by the velocities Ua = clVLa and Ub = bfls of the inner 
and outer cylinder surfaces where Qa and & b are the cor- 
responding angular velocities. As long as the two angular 
velocities are constant along the axis of the TC-system, 
a simple closed form analytical solution can be found as 
a sum of a solid-body rotation Vg ~ r 1 and a potential 
vortex Vg ~ r~ 1 ; see for instance Ref. [ID] . 

In this work we are concerned with the interaction zone 
between two co-current Couette flows which arises if the 
annular velocity f2 of one or both of the cylinders changes 
abruptly as a function of z from tt~ to fl + at z — 0. The 
velocity of the inner cylinder, for instance, varies discon- 
tinuously from U7 = aSl^ at z < to U\ — aVl\ for 
z > 0. Sufficiently far away from the discontinuity we an- 
ticipate that the azimuthal velocity will depend only on 
the radial coordinate r. In the vicinity of z = 0, on the 
other hand, Vg — Vg(r,z), and the solution is governed 
by the Laplace equation ([!]). In order to solve this ellip- 
tic equation, boundary conditions are required not only at 
the cylinder surfaces but also in the axial direction. To fa- 
cilitate the derivation of an analytical solution we impose 
axial periodicity, i.e., 

Vg(r,z) = Vg(r,z + 2L) (3) 

where 2L is the periodicity. At the end of the analysis 
the limiting case L — > oo will be considered such that the 
solution satisfies Neumann boundary conditions as z — > 

±00. 

3. Analytical solution 

By the standard method of separation of variables, as- 
suming 

V e (r, z) = R(r)Z{z) 




Figure 1: Top: Asimuthal velocity Vg as function of r and z for 
b = 2a = 1, = Ug = 1 and [Tt = Ug = —0.5 (arbitrary units). 
Bottom: Absolute difference between the flow in the top panel and 
the reference flow, Eq. plotted as equidistant shaded bands to 
illustrate mixing layer shapes with different SU (sec section [3. 2|. 



gives 

z , x _ f M + B z, n = 

1 A n cos k n z + B n sin k n z, n>l 

RM _t CV + A^, n = 

n[r) \ C n h{k n r) + D^Kr), n>l 



where 

T17T 

k„ = jr-, neN (4) 

and An, B n , C n , D n are undetermined constants and No is 
the set of positive integers, including zero. 

We may use standard theory of Fourier series to deter- 
mine the coefficients A n and B n to satisfy the boundary 
conditions along the split cylinder boundaries. Combined 
with the remaining boundary conditions one readily ob- 
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tains the analytical solution as a Fourier series: 



V e ={U A ) 



b/r — r/b 



r/a — a/r 



b/a — a/b ' s Dl b/a — a/b 

K x {k n r)h{k n b) - K^KtyhiKr) 



, f AU K^k n r)h{k n l 
£j\ A K x {k n a)h{k n l 



AU 



% b) - K x [k n b)h{k n a) 
Ki(k n a)h(k n r) - ifi(fc„r)7i(fc„a) 



Ki(k n a)Ii(k n b) - ifi(fc n 6)/i(fc„a) 
nn 



sin(fc n z) 



where 



AU A , B =U+ B - U a ,b- 



(5) 

(6) 
(7) 



Ii and Ki are modified Bessel functions of order 1. One 
sees immediately that the first two terms (which come from 
the n = term) correspond to the case of normal TC flow 
(without cylinder splitting) in which the inner cylinder 
rotates at angular frequency (Ua)/ol and the outer with 
(Ub) /b. The terms proportional to AUa.b then contain 
the breaking of translational symmetry along the z axis, 
and we refer to these as the asymmetry terms. 

In the special cases where either just the inner cylinder 
or just the outer cylinder is split, we get velocity profiles 
depending on a single parameter, generalising that consid- 
ered in Ref. [5] , respectively, 



f] [1 - sin(fc ra z) 

* J 11 



V e = b/r-r/b | K A 
cl(£Ia) b/a — a/b tt 

v ' ' ' n— 1 

x gi(fe»r)/ 1 (fc n b) - ^ 1 (fc n 6)J 1 (fc ra r) 
Ki(k n a)Ii{k n b) - K 1 (k n b)I 1 (k n a) ' 



(8a) 



and 



Ve r/a — a/r 



b(n B ) b/a- a/b 



7T 



E 

n=l 



[1 -(-!)»] 



sin(fc„z) 



K 1 (k n a)I 1 (k n r) - K 1 (k n r)I 1 (k n a) 
K 1 (k n a)h(k n b) - Ki(k n b)h(k n a) ' 



(8b) 



where now (CIa,b) = + anci ^a,b = (^a,b~ 

b)I (^a.b)- The full solution ([5| may be seen to simply 
be the sum of (8a) and (8b), which in fact follows from the 



linearity of the governing equation. 

In the limit L — > oo the sum over n becomes an integral, 
and the case of Neumann boundary conditions at z — > ±oo 
is obtained by means of e.g. the Euler-Maclaurin formula 
(e.g. Ref. [UJ, §23.1) as 



V e =(U A ) 



b/r — r/b 



(Ue 



r/a — a/r 



AU 



b/a — a/b b/a — a/b n j 

- K x {£b)h{£a) 



d£ 



sin(£s) <! AU A 



(9) 



and the special cases of either just inner or just outer cylin- 
der split, respectively, 



and 



Ve _ b/r-r/b A f°° d£ . 
a(Q) ~ b/a -a/b + tt J £ 



Ve _ r/a-a/r A f 00 d£ . 
b(Q) ~ b/a -a/b tt J £ ' 
x K^a)h^r) - K^h^a) 



(10a) 



(10b) 



For reference in the following we will use as reference 
flow the (unphysical) velocity profile made up of two iso- 
lated creeping Taylor-Couette flows which meet discontin- 
uously at z = 0: 



ref _ b/r- r/b fTT+ ^ 

b/a-a/b [UAW ' j 1 ' 



vr = 



: [U+Q(z) + U^Q(-z)} 



r/ " " /r \u+e{z) + u^e(-z)] (iia) 



b/a — a/b 
b/r — r/b 



b/a — a/b 
r/a — a/r 



b/a — a/b 



(U A ) + iSg(z)AU A ] 
(U B ) + tSg(z)AU B ] 



(lib) 



where &(z) is the unit step function and Sg(z) is the 
signum function. The latter form is particularly useful 
in the following for comparison with Eq. 

An example flow is shown in figure [T] We observe how a 
mixing layer forms between the two Taylor-Couette flows 
whose width vanishes at either wall discontinuity. 

3.1. The limit of plane co-current Couette flows 

As a test of our solutions, the expressions obtained in 
Ref. [5] for the case of a mixing layer between two plane 
Couette flows should be regained in the limit where 

The large sum in Eq. ([5| receives its main contributions 
from n such that k n < 1/(6 — a), since the large frac- 
tion decays exponentially for larger values of k n . Then 
k n a, k n r, k n b 3> 1 and we may use (e.g. Ref. [TT] §9.7) 



H*) 



12-KZ 

so that, with rj = r — a, 



Ki(z) 



Ki(k n a)h(k n r) 
K 1 (k n r)I 1 (k n a) 



2fc nV / ar' 

e k n r) 

2k ny /ar J 
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etc., to yield 

ifi(fc„g)Ii(fc n r) - K 1 (k n r)I 1 (k n a) 
Ki(k n a)h(k n b) - Ki{k n b)h{k n a) 
K 1 (fc»r)/ 1 (fc»b) - Kxjkr^hM 
Ki{k n a)Ix{k n b) - Kx[k n b)h{k n a) 

where we have used that \fbjr ~ 
additionally that 



sinh k n r\ 
sinh2£;„ft' 
sinh k n (2h — r/) 
sinh 2k n h 

'rja ~ 1. Noting 



r/a — a/r r\ b/r — r/b 2h — T) 



b/a — a/b 2h b/a — a/b 



2h 



we obtain the planar result with periodic boundary condi- 
tions in the spanwise direction 



V e -*U 



2ft 



2h 

OO 



'-(Ua 



2ft 



(Ue 



2AU A ^ sinh 



(2n+l)ir(2h-z) , (2n+l)ny 
L bU1 L 



2AU f 



sinh 



(2n+l)irz 



2n + l 

{2n+l)-Ky 



n=0 ° Um L 



2n+ 1 



(13) 



where, for comparison with Ref. [S], we define y = z and 
2h = b — a and z = r — a. In the special case (J/a) = 
A {/a = this is exactly the result reported in Ref. [5]. 

In the case L — > oo one obtains the analogous Neumann 
boundary condition result [9] 



U(y,z) = (U A ) 



2ft 



2h 



(Ue 



2h 



AU A 



+ 



AC/, 



B 



arctan 



arctan 



t.anh I ~1 tan ( % '"" 1 



7T 



4ft 



, . ivy \ I uz 
tanh — - tan — - 
4ft 



(14) 



Mixing layer width 
An approximate analytical expression for the spanwise 
width of the mixing-layer will now be derived, similar to 
the analysis in Ref. [3]. Let us define the mixing-layer so 
that it covers the region within which the flow Vg{r, z) in 
Eq. ([9]) differs by more than some small measure SU from 
the discontinuous (and unphysical) flow profile in Eq. (Ill 



made up of the two isolated Taylor-Couette flows jumping 
sharply at z = 0. For numerical purposes we choose SU = 
<5Max(|A£/4|, |AZ7b|) with 5 a small parameter which we 
take to be 0.01 in the numerical examples. The edges of the 
mixing layer, which we call z m { x , then solve the equation 



SU 



I f{-"} Sin ^ 



AUAjK-n-^ + AUB^-rj- 1 ) 
A -A- 1 



(15) 



where Sg is the signum function, we make the shorthand 
definitions 

b r b , . 

A=-, i]=-, n=- (16) 



and {•••} is the expression in curly braces in Eq. Q. 
One may show that the mixing layer is always symmet- 
rical about z = 0. 

Noting that the integrand is an even function of £ we 
may expand the integration range to cover the whole real £- 
axis and close the contour in a large semicircle in the lower 
half-plane, noting that sin£z = — Sg(z)Im{exp(— i£|z|)}: 



- / -H---}sin£z 

T JO £ 

_Sg(^) x „ /d£ 



.Mf) Im r f { ... }e -^i 



2tt 



Tm 



{•••} 



=Sg(z)ReV Res-f.-le^l 

3=0 



(17) 



where T is the integration path previously described, j 
sums over the poles £j of the integrand, and the prime 
on the summation mark signifies that the pole at £ = 
is counted with half weight. These poles lie at £ = 
and along the imaginary £-axis and the ones along the 
negative imaginary axis we denote £j = — iyj. Calculating 
the residue at £0 = we find that it exactly cancels the 



second term on the right-hand side of Eq (151 



For the remainder of the sum we make the approxi- 
mation of keeping only the first term, corresponding to 
£1 = ayi which is the smallest of the values £„ that solve 



Yi(Q Ji(Q 



(18) 



where J n and Y n are Bessel functions. We then obtain the 
solution 



where 



no 



a , SU 



2 au aPi ( v C;K) + ^u BPl (c,vO 



AC 2 p 2 (C,AC)-po(C,AC) 

Pn{xi,X 2 ) =Jn(xi)Y n (x 2 ) - J n {x2)Y n (xi), 

and we have used (Ref. [H] §9.1) 
d 



(19) 

(20) 
(21) 



d.r 



Pi(ax, bx) = —^abx[p2(ax,bx) — po(ax, bx)]. 



We can obtain a better approximation by including one 
further term of the infinite sum of residues and use that 
£2 ~ 2£i [for large £ the solutions to (18 1 are asymptoti- 



cally £ n ~ mr/(X — 1)], which we will refer to as a second 
order approximation: 



, , a , g i^ 2 (Ci)+4^ Sl F(£ 2 T-^(£i) 
l^mixl- Ci iog 2 ^ (C2) \ZZ) 

with si = Sg[F(£i)]. This tends to the first order expres- 
sion (fl9| when SU — > 0. 



a) z/a Only inner cylinder split 



b) z/a Only outer cylinder split 
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Figure 2: Mixing layer for the three cases of only inner cylinder split, only outer cylinder split, and both cylinders split. The non-split cylinder 
is at rest (but this makes no difference to the mixing layer width). Velocities in arbitrary units. The exact boundary layer is the shaded 
area and the dashed line shows the first-order analytical approximation, Eq. (19} . In panels c, d and e also the second order approximation 
Eq. i|22[) is shown as a dotted line. 6U = 0.01 max J 4,s(AC/) and b = 2a = 1 in all panels. 



Mixing layers for a few different cases are shown in 
Fig [2] The shaded area represents the mixing zone cal- 
culated from the exact analytical solution in Eq. ([5]), or 
if appropriate, from Eq. Q. Outside of the circumfer- 
ence of the shaded area, the exact solution deviates from 



the approximate discontinuous solution in Eq. (Ill by less 



than 1 per cent. When both cylinders are split, two-lobed 
mixing layer shapes are observed, as shown in panels c-c 
in the figure. In this case the second order approxima- 



tion (22 1 follows the shape of the greater of the lobes with 



excellent precision whereas the first order approximation 



(19 1, although taking the shape of a single lobe, gives a 



reasonably good outline of the contour of the mixing layer 
albeit with somewhat limited accuracy when the two lobes 
are of similar size, as exemplified in panel d. The maxi- 
mum width of the mixing layer is however well predicted 



by Eq. (19) 



4. Shear stress and torque on rotating cylinders 

It is of interest to investigate how the mixing layer 
changes the viscous torque experienced by the rotating 
cylinders as compared to the standard Taylor-Couette set- 
up in the creeping flow regime. The shear stress acting on 
the cylinder surfaces at r — a, b directed along the local 
6*-axes are obtained via the rates of strain from Newton's 



law of viscous friction, 



t a = ixd r V e 



t b = -fJ.d r V e 



(23) 



with fi being the dynamic viscosity. The torques acting 
on the walls in the standard Taylor-Couette set-up are 
easily calculated, and generalizing slightly to the reference 



flow in Eq. (11 1, the z component (being the only nonzero 



component) is 



M 



±,ref 



±,ref 



B 



-1-K\iaL 



:2irfibL 



2U% - (A + A- 1 )^; 

A - A- 1 
2U% - (A + A- 1 )^ 

A- A- 1 



(24a) 
(24b) 



Here 2L is the total length of the split cylinders. 

It turns out that the additional torque on the individual 
halves of a split cylinder due to the presence of the mix- 
ing layer is formally infinite. Mathematically this can be 
seen from the fact that d r Vg at r = a and b has a simple 
pole-like singularity ~ 1/z at z = 0. An integral over the 
cylindrical surface area integrating to z = must thus nec- 
essarily be logarithmically divergent. The same is found 
to be true of the shear force in the planar case considered 



in Ref. [9], Eq. (14). The radial derivative of Vg is shown 
m Fig. [3] for the same parameters as in Fig. [l] Clearly, this 
infinity is not physical, a point to which we shall return. 
We will discuss the behaviour of d r Ve and the local shear 
stress further in the following. 
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0.2 0.0 0.2 0.4 0.6 



Figure 3: Radial derivative of Ve with same parameter values as in 
figurc^and the contours of the same (bottom). The radial derivative 
has a pole singularity at the two discontinuity points. 



An explicit expression for d r Vg is trivial to write down 
from Eq. ^ . We are interested in the difference in viscous 
shear resistance due to the mixing layer, compared to that 



of the reference flow field (11), as given by Eqs. (24), and 



subtracting the derivative corresponding to the reference 
case the difference may be written 



Ad r Vg 



1 ) -AU a {k + k 



Sg(» rAU B (r) + V~ 



2r 



A-A- 



sin£z 

(25) 



It is convenient in the following to move the last term into 
the integral over £ by noting the Fourier identity 



which is valid for any real z. 



d£ 



sin£z 



At the inner and outer cylinders Eq. (25) reads, explic 



itly 
Ad r V B 



7T J { £a(A-A ) 

AU B /ta + AUAlKii&hitb) - K^b) ' 



(26a) 



and 

- Ad r V 



d£ sin £ z 



1 

r=b 7T jq 

AUA/tb + AUBlKiiZtyhfca 



AU B (\ + \- 1 )-2AU A 

^(A-A- 1 ) 
-*i(£a)Jj(#) 



K^hitb) - K x {ib)h{ia) 



(26b) 



We may use the path integral technique as in the previ- 
ous section to write this as a sum of residues. The ma- 



nipulation (17) holds also for the integrals in (25) and 



(26). Calculating the residues, which are found at the 



same poles as in the previous section, we again find that 
the residue of the pole at z = is exactly cancelled be- 



tween the two terms in the integrands of Eqs. (26a) and 



( 26b ) . The derivatives may then be written as sums over 



the remaining residues, 



Ad r V t 



-CA z \/ a 



CfMO.ACO-PoCC-.ACi)] 



(27a) 



and 



-Ad r Vs 



r—b 



<j\ z \/ a 



with 



aSg(*)xy 

x au b cm^iKj)~^Ua/\ 

Cj [P2 (Cj , AC, )-Po{(j, Ki )} 



Qn{x,y) = J n (x)Y^(y) - j' n (y)Y n (x). 



(27b) 



(28) 



The values Q are still solutions of Eq. (18). When \z\ is 



of order a or higher, these sums converge rapidly and are 
numerically much cheaper than evaluating the integrals. 

As already mentioned, the torque contributions to the 
half-cylinders from the singular points are formally infi- 
nite when the shear stress is integrated over the cylinder 
surfaces all the way to z = 0. Let us therefore exclude a 
thin band of width Az close to the discontinuity and cal- 
culate the torques on the upper (+) and lower (-) cylinders 
according to 



M\ =2TT^a 2 / dzd r Ve 



Az 



-Az 



dzd r Vg 
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and similarly for the outer cylinders. 



With expressions (27 1 for the strain rates the corre- 
sponding torques on the upper (+) and lower (-) inner 
(A) and outer (B) cylinders then take the form 



AM 



(<>) 



± 



Z \ -\ LSj 



AUb ] 
AU A i 



-Qj Az/a 



(29a) 



and 



AM 



AM 



(0) 



± 



2 ~ fcq^XQ-^e-^*'* 
X U C^ 2 (G,AO)-po(0,AO)] ' 



We have defined the reference torques 

AM^ 0) = 2n^a 2 AU A , AM^ = 2n^b 2 AU E 



(29b) 



(30) 



As Az — > the sum is dominated by large values of 
Q rj 7rj/(A — 1), for which, using asymptotic expansions 
of Bessel functions (Ref. (TT] §9.2) 



2 C,f/i(Cj,Ag) 



AU B 
AUa 



A C?[P2-Po](G,AG) 
2[C J 9i(0,A0)-^ 1 



AA£/ B 



Po](Cj,AO) 
where we noted 

?i(G,ACj) ~gi(AG,G) 



A AC/b 
2 1 J Al^ 



1 



2a/A 



(-1> 



AC/r 



tt^VA 



cos[Cj(A — 1)] 



[P2 - Po](Q,Kj) ~ - cos[G(A - 1)] 

for large Q. We have inserted (j ~ ^/(A — 1) for large j. 
The shorthand notation introduced in the denominators 
should cause no confusion. We obtain 
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(31b) 



for Az — > 0. 

In the mathematical formulation of the problem in Sec- 
tion^ the angular velocity of either the inner or the outer 
cylinder (or both) was assumed to exhibit a discontinuous 
change at z = 0. This discontinuity gave rise to a pole-like 
singularity when the wall shear stresses and the accompa- 
nying torques were calculated in Section |4j In practice, 
however, a discontinuous change of the rotation rate is a 
mathematical artifact which cannot be reproduced in a 
laboratory TC apparatus. If the inner cylinder, say, is 
divided in two differently rotating halves, a tiny gap will 
exist between the two parts and tend to smear out what 
was meant as a discontinuity in the angular velocity of the 
cylinder. 



5. Concluding remarks 

An analytic solution has been derived for the shear layer 
between two co-current Taylor-Couette flows. A steady- 
state solution which obeyed axial periodicity was obtained 
first and thereafter simplified to an infinitely long period. 
In the limit of large radii and a narrow gap, the earlier 
solution by Narasimhamurhty et al. [9] for the bilateral 
shear layer between two parallel plane Couette flows was 
recovered. 

While co-current plane Couette flows are challenging to 
realize in a laboratory experiment, a Taylor-Couette flow 
apparatus can be constructed such as to allow the inner 
cylinder to rotate with an almost discontinuous angular 
velocity. This set-up will enable investigations of, for in- 
stance, particle dispersion or chemical reactions in a bi- 
lateral shear flow with shearing motion not only in the 
(r, #)-plane both also in the perpendicular (#,z)-plane. 
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